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I. INTRODUCTION 



A successful phenomenology of the hadron spectrum has been obtained using non- 
relativistic potential models, which tentatively simulate the low-energy limit of QCD Jl] [| . 
The interquark potential usually contains a linear part which describes QCD confinement 
and is supplemented by a Coulomb term which may be attributed to one-gluon exchange. 
Spin-spin, spin-orbit and tensor terms are added, analogous to the Fermi-Breit components 
of QED potential, derived from v/c expansion. 

There are obvious difficulties. Large relativistic corrections can be anticipated for light 
quarks. The naive superposition of a Coulomb and a linear term may be too schematic. 
Results rely on the potential at intermediate distances (0.1 fm < r < 0.5 fm) where neither 
the perturbative nor the string limit holds. Nevertheless, the success of potential models 
indicates that somehow delicate relativistic and field-theory effects are hidden in the pa- 
rameters. As none of the more ambitious approaches, for instance lattice calculation Q, is 
yet able to produce very precise results, it is still justified to use potential models as tools 
to analyze hadron properties, with the hope of better understanding the non-perturbative 
limit of QCD. 

Among the observables of interest, isospin- violating mass differences have retained much 
attention. Earlier studies of these mass differences || have been reconsidered within 
constituent-quark models. In general, the n — p, S~ — E°, S~ — E + , H~ — S° splittings 
of the nucleon, £ and S multiplets are well reproduced, this fixing the quark-mass differ- 
ence Am = ma — m u . Predictions for charmed baryons can then be supplied. Some results 
concerning the S c and S c multiplets are shown in Table [I], together with experimental data 
||. The estimates of Wright @ and of Deshpande et al. || are not really potential models, 
they are shown only for information and comparison with the others |9|-|T7| 

Some of the models include only a fraction of the possible contributions, for instance 



2 



electrostatic interaction is accounted for, but the mass dependence of the chromomagnetic 
interaction is neglected when replacing a d quark by a u one. This is hardly justified. As 
underlined, e.g., by Isgur |16|], these isospin splittings arise from several canceling contri- 
butions, so that each effect should be carefully computed and even small terms should be 
incorporated. 

The most striking feature of Table | is the wide spread of predictions. Next come the 
observation that none of the models is compatible with the presently available data 0. In 
particular the predicted — splitting tends to be smaller than the PDG average 6.3±2.1 
MeV |J. However, a preliminary result by CLEO indicates a smaller value S° — = 
2.5 ± 1.7 ± 1.1 MeV |T8|] . This collaboration has also detected candidates for the internal 
spin excitations and E'^j" (total spin 1/2 with the light-quark pair sd or su mostly in a 
spin-triplet state). Their measurements indicate a splitting S'° — E7^~ ~ 1.7 MeV with a 



large error bar. There are also data from CLEO on the 5* states [[18 



The S c multiplet is the most puzzling. The S+ + — TP C splitting is usually larger than 
— while data seemingly favour the reverse. In other words, most models predict an 

ordering of and E° which is not seen, to the extent one can draw any conclusion 

from the data. 

The present investigation is motivated by the discrepancy between data and models. We 
wish to understand whether this problem points out a general limitation of potential models, 
in particular those based on one gluon-exchange, or can be solved by reconsidering the choice 
of parameters and removing unjustified approximations in the three-body problem. 

For this purpose, first we carefully estimate the role of each contribution to the splittings 
within specific potential models. This should measure to which extent previous calculations 
suffer from neglecting some effects or treating them approximately. Then we analyse the 
sensitivity to the choice of potential, to see whether fitting the data can be achieved by 
an appropriate tuning of parameters or it is out of reach of this approach. Predictions are 



listed for a number of isospin multiplets, making possible a comparison with all available 
experimental data, and with those one could expect to be measured in the near future. 



II. MODEL CALCULATION 

A representative quark model is proposed in Ref. ||, where the potential is tentatively 
designed to fit both meson and baryon spectra, using the empirical rule 

V QQQ (r 1 ,r 2 ,r 3 ) = ~Y / V(\? i -? j \). (1) 

i<j 

In 0, the quark-ant iquark potential reads 

V(r) = — - + Xr p - A + ~6(r, r ) ff x ■ a 2 , (2) 

with 5(r, ro) = exp(— r/ro)/(47rrQr) being a smeared form of the contact term. The param- 
eters are (in the units used by the authors [fj]): 

p=l, K = 102.67 MeV.fm, «' = 6 x 102.67 MeV.fm, 
(BCN) A = l/(0.0326) 2 MeV/fm, A = 913.5MeV, r^ 1 = 2.2fm" 1 , (3) 
m q = 337, m s = 600, m c = 1870, m 6 = 5259 MeV. 

One of the difficulties in the above model is that the spatial extension ro of the spin- 
spin term is too large to describe a short-range interaction between heavy quarks. As a 
consequence, the J/ty — rj c hyperfine splitting is not well reproduced. This is why Ref. ||, 
following for instance Ref. |[L9|| , introduces a flavour dependence in ro, namely 



(2/72 Tfl ■ \ 
— H" , (4) 
rrii + in j J 

while a Gaussian form 5(r, ro) = exp(— r 2 /rQ)/(7r 3//2 rQ) is now adopted. Both models AL1 
and API of Ref. fit very well the meson spectrum. The parameters are 



p = 1, K = 0.5069, k' = 1.8609 A = 0.1653 GeV 2 , 
(AL1) A = 0.8321 GeV, B = 0.2204, A = 1.6553 GeV s_1 , (5) 
m q = 0.315, m s = 0.577, m c = 1.836, m fc = 5.227 GeV, 



and 



p = 2/3, k = 0.4242, k' = 1.8025, A = 0.3898 GeV 5/3 , 
(API) A =1.1313 GeV, 5 = 0.3263, A = 1.5296 GeV 5-1 , (6) 
m g = 0.277, m s = 0.553, m c = 1.819, m b = 5.206 GeV. 

None of these models include tensor forces, since this interaction is not expected to give 
important contributions in hadron spectroscopy, at least for ground states . 



For studying isospin breaking, we have allowed for m u ^ nid and added to the potential 
the electrostatic interaction between quarks. The difference Am = m^ — m u between d and u 
quark masses has been adjusted to reproduce neutron-proton and S _ — S + mass splittings. 
In particular, we have taken m u = 327 MeV and = 338 MeV for AL1, m u = 337 MeV 
and m d = 353.85 MeV for BCN and m u = 277 MeV and m d = 300.5 MeV for API. Our 
estimate of Am for constituent quarks is larger than the common wisdom for current quarks, 
Am ~ 4 MeV. However, it must be noticed that dressing quarks modifies this quantity, for 



the cloud of virtual states depends on the flavour of the quark it is surrounding Plj^2 . 



Baryon masses have been obtained using two reliable numerical methods p3| , |24"|| . In the 
first case every contribution was included non-perturbatively in the variational procedure, 
while in the second one the electromagnetic terms were treated perturbatively. The perfect 
agreement of two results indicates that a good convergence has been reached. 

The various contributions to n — p and to charmed baryons mass differences are shown 
in Table [Q|, for the specific model ALL As hinted previously delicate cancellations occur, 
requiring an accurate treatment of each term. In a flavour-independent potential the energy 



of a given state decreases when any of the constituent masses is increased, i.e. (AT + AW + 
AB) I 'Am < in the notations of Table [H]. This is observed in our calculations, though in 
the S c case, the flavour dependence of spin-spin term goes in the opposite direction. As for 
the kinetic energy itself, one expects [[26|] AT/ Am > if confinement dominates the binding 
process, and AT/ Am < when the Coulomb part becomes more influential. 

Before discussing the results obtained for baryons, we have investigated the splittings 
among mesons, for instance D + (cd) — D°(cu). The results are shown in Table |T1|, for the 
three models BCN, AL1 and API. Am is adopted to reproduce m n —m p in the baryon sector. 
An acceptable agreement is found for AL1, while BCN and API seem to be disfavoured, 
leading to a larger overestimation of the K° — K + splitting. Actually, AL1 is somehow an 
improvement of BCN, while the confining part oc r 2 / 3 of API does not agree with lattice 
results and fails for heavy mesons as well. In the following AL1 will be our benchmark, the 
results for API and BCN will also be given in order to show how stable are the predictions 
respect to the choice of the potential. 



Our results for baryons splittings are shown in Table IV 



A check of the consistency of our calculation is provided by the relations listed by Franklin 



25fl . When one switches off the electromagnetic interaction, one deals with energies which 
evolve continuously when going, for instance, from (Quu) to (Qdd) via (Qud). Hence a mass 
combination 

8{E Q ) = (Quu) + (Qdd) - 2{Qud) (7) 

receives contributions mostly from Coulomb effects. If the latter are treated at first or- 
der with a wave function (Qqq) properly averaging that of (Quu), (Qdd) and (Qdu), then 
<5(£q) ~ a(r~q), i.e. the charge of heavy quark Q disappears. 

If furthermore the (qq) part of the of the (Qqq) wave function does not depend much on 
the mass of the heavy quark Q and on the coupling of the spin of Q to the (qq) spin triplet, 
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then 6(Eq) should be approximately the same for S, E*, S c or multiplets. This is again 
rather well verified in our calculation. 



The possible dependence of the (qq) distribution on the mass of the third quark Q is 
discussed by Rosner [|(J as a "three-body effect" . It was investigated previously f27|j28| in 
the literature, not for (r~g), but for the somewhat similar matrix element (5^(i qq )) that 
enters the calculation of the hyperfine splittings in usual quark models. In an approach a la 
Breit-Fermi, the ratio 

2S- + S-3A 

R ~ 2A-2N (8) 
reveals the ratio of the (qq) short-range correlations in (sqq) and (qqq). Similarly the ratio 

R ' = §r| P) 

gives a comparison of (qs) correlations in (ssq) and (sqq). The experimental values R ~ 
1.04 and R' ~ 1.12, as well as the detailed three-body calculations p7| , p8| show that, as 



conjectured by Rosner [2^|, the q\ and q2 quarks tend to bind more intimately within (qiq^qa) 
when g 3 becomes heavier. The effect is about 5 to 10 % when q 3 changes from ordinary to 
strange, may be slightly more from strange to charmed. In would be desirable to reach a 
deeper understanding of this property, beyond numerical investigations. A possible starting 
point is given by the harmonic oscillator 

H = v\ + vl + ap 2 3 + r\ 2 + (3(r 2 13 + r 2 2 3 ) , (10) 

where the r 12 dependence of the wave function factorizes out and is easily shown to be 
independent of the inverse mass a (but it does depend on the strength ft). 
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III. DISCUSSION 



The isospin-violating splittings of light and heavy baryons are shown in Table [IV]. The 
S~ — S° and A — A ++ splittings and even (within large errors) the electromagnetic split- 
tings for the excited states S*(1385) and S*(1530) come out in good agreement with the 
experimental data ||. Some problems appear however for charmed baryons. 

More precisely, while the experimental datum X+ + — S° = 0.8 ± 0.4 MeV is well repro- 
duced, one finds an, albeit small (of the order of —0.5 MeV), negative E+ — E°, at variance 
with the experimental datum £+ — S° = 1.4 ± 0.6 MeV. Also the result for — S° = 2.2 
MeV is smaller than the, still rather imprecise, experimental datum — 5° = 6.3 ± 2.3 or 
4.7 ±2.1 MeV, where the first number corresponds to the particle data group average and 
the second to their fit. The problem is not solved using API or BCN. 

To summarize at this stage, the splittings of charmed baryons do not agree with experi- 
mental results, when they are calculated from potential models supplemented by electrostatic 
forces and a mass difference Am between d and u quarks. 

Reasonable changes of light quark masses do not modify substantially this situation. 

An effect which we have neglected up to now is the electromagnetic dipole-dipole interac- 
tion between quarks, whose dominant term (neglecting the small contribution of components 
of the wave function with non-vanishing angular momentum) is: 

2 7i qiqja 



3 mirrij 



S^i^Bi-Bj (11) 



where are quark charges in units of electron charge. Albeit surely present, one expects it 
to be smaller than usual Coulomb interaction. Our numerical results including this dipole- 



dipole term are shown in IV. We used a regularized form 5 for 5(r) } taken to be the same 



as for the strong spin-spin force. The magnetic contribution goes in the right direction, 



S 



but remains too small to push the computed masses significantly closer to the experimental 
ones. 



This problem raises the question whether some contribution has been forgotten. For 
example, in some models adjusted to reproduce meson and baryon masses simultaneously 
introduce in the baryon sector an ad-hoc 3-body term of the form m 

D 3 + - 3 - h . (12) 

This parameterization is purely empirical. For the AL1 model, the parameters are D 3 = 
0.07376, A 3 = -0.05546 and b 3 = 1/4. 

As it depends on masses, this term gives a contribution to isospin breaking effects as well. 
The 3-body term flT2| ) slightly improves the description of the electromagnetic splittings of 
light baryons. However, how is evident by inspecting Eq. (|i~2"D, the contribution — S° 
goes in the wrong direction. When the 3-body term ( |l"2"D is accounted for one obtains 
typically £+ - £° ~ -0.7 MeV (instead of -0.55 MeV), and E° c -E+ ~ 1.5 MeV (instead of 
2.58 MeV), with little dependence on the choice of parameters. Of course, one could think 
of more complicated three-body interactions, but their form remains completely arbitrary 
and somehow the appealing features of potential models are lost once one violates flavour 
independence and gives up the link between quark-quark and quark-antiquark forces. 

Another possibility which can be explored is the running of a s , which leads to a reduced 
coupling when heavy quarks appears for the scale is chosen to be proportional to the masses 
involved (of course problems related to the precise choice of the scale and to the unknown a s 
behaviour at small scales emerge). Such an effect would decrease the strenght of the spin- 
spin term involving heavy quarks, but this would not go in the right direction for changing 
the order of S c states. 

Finally, we discuss now two interaction terms which have been contemplated in addition 
or in replacement of chromomagnetism. 



The first comes from the contribution of instantons. The non-relativistic form of a 
potential mimicking 't Hooft interaction |29[ has been elaborated in Ref.s [3^,^]. The 



value of the coupling must however be fixed phenomenologically. Interesting results have 
been obtained on hadron spectroscopy with models including this instanton term replacing 
[ 3Ti| , |3"T| ] or supplementing [32] the chromomagnetic force. 



However, when the instantonic potential is considered only as a further correction to Eq. 
(0), it does not contribute substantially to S c mass splittings, for it is inversely proportional 
to the quark masses and vanishes for a quark pair with spin 1. Thus it cannot help solving 
the problem of S c splittings. It gives, anyway, a positive contribution, albeit not to be 
expected quite large (of course the precise numerical value will depend on the choice of the 
coupling), to — H°. 

The second type of interaction deals with meson exchange between quarks. There is a 
rich literature on the subject, which has recently be revisited by Glozman and collabora- 
tors (see e.g. Refs. |33| and references therein), who have adopted rather an extreme point 



of view where the chromomagnetic force is completely removed. These authors obtain a 
surprisingly good fit to light and strange baryons. In this approach, the study of electro- 
magnetic splittings of baryons is somewhat reminiscent of isospin violating effects in nuclear 
physics, where one accounts for a difference between n and 7r° masses and their couplings 
to nucleons. This remains to be studied. However, the extension of Glozman model to heavy 



baryons seems problematic, notwithstanding some initial attempts [35 



IV. OUTLOOK 



In conclusion, we find that non-relativistic potential models do not permit to reproduce 
the data on — S° and — S° mass splittings, despite the good agreement obtained 
for light baryons. Of course experimental data need further confirmation, and need to be 
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extended to beauty and double-charm sectors, to see if this discrepancy persists. 

Previously, Franklin pointed out that the experimental data on charmed baryons 
violate mass relations which are expected to hold within large class of quark models. We 
have checked that the mass obtained from an accurate solution of the three-body problem 
fulfill the Franklin relations. So if the difficulty persists, its solution should be searched 
in a intrinsic limitation of usual quark models, for instance in the need for new dynamical 



contributions, such as electromagnetic penguins . 

The present situation is somewhat a paradox. The A c and Q c and the average E c , 
H c states are reasonably described by simple potentials, i.e., one seemingly controls the 
behaviour of the ground-state baryons when an ordinary quark is replaced by a charm one. 
Meanwhile one does not understand the effect of a more modest move, when a up quark is 
changed in a down one. 
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TABLE I. Predictions of different models for charmed baryons electromagnetic mass splittings 



Model 


£ c ++ - s° 


£+ - £0 




Experiment Q 


0.8 ± 0.4 MeV 


1.4 ±0.6 MeV 


6.3 ±2.1 MeV 


Wright 


-1.4 


-2.0 


3.1 


Deshpande et al. || 


-(3-18) 


-(2.5 - 10) 


4.5 - 12 


Itoh § 


6.5 


2.4 


2.5 


Ono 


6.1 


2.2 


1.8 


Lane and Weinberg [p^] 


-6 


-4 


4 


Chan |11| 


0.4 


-0.7 


3.2 
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Lichtenberg 
Kalman and Jakimow [14] 



Capstick [15 



Isgur [16 



Richard and Taxil [17] 



II 



3.4 
-2.7 
1.4 
-2 
3 

-2 



0.8 
-2.2 
-0.2 
-1.8 
1 

-1 



1.1 
3.6 
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TABLE II. Different contributions (in MeV) to the neutron to proton and charmed baryons 
mass differences: Am is the change of constituent quark masses, AT the difference of kinetic 
energies, AW the variation of the expectation value of the Wigner term (independent on spin 
and isospin), AB the difference of Bartlett components (oc Oi<7j) and finally AC comes from the 
Coulomb electric interaction. Here we use the model ALL 



Baryons 


Am 


AT 


AW 


AB 


AC 


Total 


n — p 


11 


0.33 


-6.45 


-2.86 


-0.76 


1.24 




-11 


0.56 


7.55 


1.39 


1.14 


-0.35 


£++ - £° 


-22 


0.79 


15.78 


2.07 


4.53 


1.20 


ctO _ ^+ 


11 


-3.31 


-5.56 


1.72 


-1.01 


2.83 
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TABLE III. Comparison of the isospin-breaking splittings of mesons (in MeV) obtained from 
several potential models: Bhaduri et al. (BCN); Silvestre-Brac and Semay (AL1 with linear con- 



finement, API with a r 2 / 3 confinement). 

Splitting Exp. (Ref. f§) BCN AL1 API 

K°-K + 3.995 ±0.034 13.15 6.64 9.56 

K*°-K* + 6.7 ± 1.2 1.55 1.36 1.28 

D+-D 4.78 ±0.10 5.37 3.78 -0.33 

D*+-D*° 2.6 ± 1.8 2.44 2.74 -0.16 

B -B~ 0.35 ±0.29 -1.46 -1.29 -6.06 

B*°-B*- -2.04 -1.23 -5.26 



1 Notice that on four available measurements, two are negative || 
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TABLE IV. Comparison of the isospin-breaking splittings (in MeV) obtained from several po- 
tential models: Bhaduri et al. (BCN); Silvestre-Brac and Semay (AL1 with linear confinement, 
API with a r 2 / 3 confinement); Richard and Taxil where the hyperfine interaction is treated per- 
turbatively (RT II has a linear central potential, RT I a r 01 one). In the column AL1 + dd, the 
magnetic dipole-dipole interaction between quarks is accounted for, in addition to the electrostatic 
potential. Some redundant splittings are shown for the ease of discussion. 



Splitting 


Exp. [Ref] 




ATI 
ALL 


Arl 


HI 11 


rvl 1 


A T 1 I A A 

Aiji + ad 




1.293318± 














n — p 


0.0000009 




1.38 


1.16 


1.29 


1.2 


1.3 


1.24 


A+ A++ 






U.04 


n no 
U.Uo 


z.Uo 






n Qfi 
U.OO 


aO A++ 
/\ — l\ 


2.7 ±0.3 @ 


Q 01 

o.zl 


on 
z.zU 


e in 

6.10 






O £ A 

z.o4 


A- A++ 






o.U4 


O.o4 


1 1 QC 
ll.OO 






£ KK 
0.00 


v- vO 

Zj — Zj 


4.88 ± 0.08 | 


I 7.09 


r 1 

5.16 


6.0 r 


n 
Z 


4 


0.24 


Zj — Zj 


8.09 ±0.16 | 


I 1 1 no 

11. yo 


O.ZD 


1 n £V 


4 




O.0 ( 


Zj — Zj 


-4 to 4 | 


§ 


/I in 

4.1U 


l.OZ 


Z.00 


4 
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